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Abstract 

This paper considers properties of nonlinear waves and solitons of Korteweg-de Vries 
equation in the presence of external perturbation. For time-periodic hamiltonian perturba- 
(<~^ ■ tion the width of the stochastic layer is calculated. The conclusions about chaotic behaviour 

C5 , in long-period waves and solitons are inferred. Obtained theoretical results find experimental 

confirmation in experiments with the propagation of ion-acoustic waves in plasma. 
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5^ 1 Introduction 

! Recently a significant attention is paid to the studying of soliton equations under external per- 
turbations. For KdV equation, in particular, it was shown, that simple sine driven force together 
O ' with dissipation can lead to a chaotic dynamics, detected by Fourier spectrum, positive Liapunov 
O ■ exponents and phase plane analysis |^. In this context the reduction of a perturbed KdV equa- 
Q . tion in the form of stationary travelling waves was considered. The justification of the application 
0^ \ of low-dimensional chaos theory to such an infinite-dimensional system as KdV equation is in 
^ ' a numerical investigation, which proves that information dimension, calculated by Kaplan- Yorke 
O ' formula, and correlation dimension within Grassberger-Procaccia method, are both between 2 and 
3 for steady waves 0. On the other hand, perturbed KdV equation can be considered directly 
whithin numerical approach, or by perturbation technique 0]. Besides this series of work it is 
p i! worth to mention here the cases, when KdV equation is influenced by dissipation and instabil- 
ity, wihtout explicit spatial or temporal dependence 0. In such a form this equation describes 
current-driven ion-acoustic instability in a collision-dominated plasma. For the strongly dissipa- 
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■ tive case overall evolution demonstrates irregular behaviour, therefore leading to nonstationary 



and irregular soliton interactions. 

In this work we add hamiltonian deterministic perturbation to the KdV equation and study 
physical consequences of this. The reduction in the form of travelling waves is made in order 
to turn the system into a second-order ordinary differential equation. The existence of solutions 
to such an equation was studied in 0, while the construction of periodic orbits can be found 
in P] for nonresonant and primary resonant cases, and in for secondary resonances. We'll, 
however, perform the analysis of chaotic properties in this situation. It is known, that for any 
hamiltonian perturbation, which makes the system near-integrable, a stochastic layer around a 
separatrix appears. This layer is bounded by unbroken KAM-surfaces, and whithin it the system 
evolves with mixing. For dissipative perturbations, however, these KAM-surfaces are broken and 
special methods like Melnikov one should be used to determine the conditions for chaos to appear. 
So, we calculate the width of a stochastic layer, which contains long-period waves and solitons. It 
is done on the basis of Chirikov criterion for the overlap of resonances, which is widely used for 
investigation of chaotic properties in hamiltonian systems ( references on original papers in this 
field can be found in [|^). Results we obtain are as follows: solitons and nonlinear waves prove to 
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be chaotic in the meaning that in a small distance from the peak of solitons and long-period waves 
there must be a region of chaotic dynamics, where these waves acquire small irregular deviations. 
This conclusion is confirmed in experiments with ion-acoustic waves in plasma, where a weak 
splash before soliton peak was registered p. 

The outline of this paper is as follows. In next Section KdV equation is reduced to ODE 
and the unperturbed solutions of the latter are considered. In Section 3 we introduce canonically- 
conjugated variables action - angle and obtain general expressions for the criterion of stochasticity. 
After that in Section 4 this criterion is applied to KdV waves directly. Section 5 contains the 
conclusions and summary. 



2 Stationary waves 

Let's consider the perturbed KdV equation, taken in the following form: 

Ut + U-U^ + f3U,,, = V{U, Ut, t/., x-vt), (1) 

where x and t denote, respectively, a one- dimensional space coordinate and time; U{x, t) is sup- 
posed to be differentiable with respect to x and t sufficient number of times; V"(f/, Ut^ Ux, x — vt) 
is a small external perturbation. As far as we restrict ourselves on the case of hamiltonian per- 
turbation, it means that within the class of functions V we'll consider only those, which will not 
contain derivatives after integration over (x — vt) = S,- It's easy to check that the general form of 
such perturbations can be represented as: 

ViU, Ut, Ux, = fi{U, ■ Ut{x, t) + MU, ■ Ux{x, t) + f,{U, 0, (2) 

with the condition [^(f/, - c ■ /i(f/, 01 = fsiU, Ou- 

Searching for a solution of equation (|l|) in the form of a nonlinear stationary wave U{x,t) = 
/((^), we obtain the following ODE after one integration over ^ ( prime means differentiation with 
respect to ^ ): 

(3f" =vf-^ + F{fj\0, (3) 

where F denotes the primitive of the function V after stationary wave substitution. Integration 
constant is absent here because it can be turned into zero by an appropriate choice of variables. 
For convenience let's rewrite equation as the following dynamical system: 

(3x = vx — — + F{x,x,t). (4) 
Equation @ without perturbation ( F{x, i, t) = ) is analogous to the equation of motion: 

Px = vx — — (5) 
of a point with mass /3 in a potential field 

2 

U{x) = y - (6) 
Full mechanical energy of such a nonlinear oscillator is equal to 
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Figure 1: (Left) Phase plane of the unperturbed system @. (Right) Plot of the potential energy 
(|^). Here v = 1 and [3 = 1. 
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where p = (3x is the momentum of the point. Phase plane of this oscillator together with the plot 
of the potential energy is presented on the Fig. 1 . 

If if = — 1^^, what corresponds to the bottom of the potential pit, then the phase curve shrinks 
to a point. In the case — < if < we have closed phase curves representing cnoidal waves of 
the form: 

x'^it) = a + h- cn'^{kt\m), (8) 

where cn{z \ m) denotes the Jacobi elliptic cosine function of modulus m. Substitution of (§) in 
d) gives: 

x™(t) = V + v[{l — 2m) + 3m ■ cn'^{kt \ m)] / V — m + 1 
T'^ = 4:K{m)/k (9) 

k = \^v/ (5{m^ — m + 1)~3 

Here K{m) is the complete elliptic integral of the first kind, T™ is the period of the corresponding 
solution, and m is the modulus of elliptic functions. Under H = {m = 1), what corresponds to 
a separatrix on a phase plane, solution (|^) yields the homoclinic orbit: 

cosh^ (±) 

where A = 2^ [|. 

3 A criterion for chaotic motion 

Let's transform variables {x,p) to the canonical pair action - angle (J - action, 6 - angle) using 
standard formulae: 

I = — Ip{x, H)dx = 1(H), e = s(x, I)= f p(x, H)dx, (11) 

2tt J ol J 

where S{x, I) is the reduced action. The Hamiltonian of a perturbed motion we can write in the 
form: 

H{i,e,t) = Ho{i) + vii,e,t), (12) 



where Hq{I) is the Hamiltonian of the nonperturbed system from and V{I,6,t) corresponds 
to the perturbation term in equation (^. 

Assuming perturbation V to be periodical on time with frequency z/, and taking into consider- 
ation the fact that nonperturbed motion is integrable, we may expand the perturbation into the 
following Fourier series: 

vii,e,t) = lEVkiii) exp[t{ke - lut], 

k,i (13) 

Vki = 

The perturbed equations of motion in variables (/, 6) take now the form: 

/ = -ii ^ kVkiil) exp[i{k9 - M] 

(14) 

e = uj{i) + ltE^exp[t{ke-iut)] ^ ' 
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where a;(/) = is the frequency of the nonperturbed motion. Let's consider the motion in the 
vicinity of one particular resonance: 

muj{Irnn) — nu = 0. (15) 



The width of this resonance on frequency is = y4Vo|ti;'|, where Vq = \Vmn{Imn) and u' = 
duj{Io)/dI. The distance between adjacent resonances is equal to 

^^al3 = \^^{Im±a,n±f3) — t^(-^m,n)| , (16) 

where a = 0,1; (3 = 0,1. As a condition of chaotic motion, which means an appearance of mixing, 
we'll use the Chirikov criterion for the overlap of resonances [|To|] : 

where Au is the minimal possible value of AcUq,^ from ([T6|). For the cases {a, (3) = (1,0), (0, 1) 
and (1, 1) calculations similar to give respectively the following conditions for the border of 
stochasticity: 
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4 Chaos of Korteweg-de Vries waves 

To apply these results to the system, governed by KdV equation, one must first evaluate the exact 
expressions for uj and u) , which appear in the conditions (|l^). The frequency u) is equal to 

27r -nk 

For convenience let's rewrite it as: 



a; 



K{z) \ 2pV3 



c sm f , ^ 

^ (20) 
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where we introduced by the correlation cos | = 1 + 3^, and z = | — ^ cot |. To evaluate J , 
one should represent it as 

/ du) dH duj ( dH\ ^ 

=dHlT = ''d^[-d^) ■ ^^'^ 

Tedious calculations give finally 



UJ 



16c^/^ sin <^K^{z)\ 2/3 sin 
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-4 cos |f (1.5; 0.5; 1; z) + 



+ fy3sin| + cos|) F(1.5;1.5;2;2) 



(22) 



Substituting (pOD and (p2D in ([T8|) one makes sure that with approach to separatrix K ^ oo under 
arbitrary small Vq. This means that for any external perturbation, however small it is, Chirikov 
criterion is executed starting with some energies, and this leads to the formation of a corresponding 
stochastic layer. The transition from regular to chaotic motioncan be found approximately from 
the condition K ^ 1. This border is defined on energy by the inequality: 
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-f^min < -f^ < 0, i^min = y (cOSV^min - !)• (23) 



Here the value of v^min is found from: 

4f ^ sin ipmin tan(v9min/3)-ft'(^mm) 



-4F (1.5; 0.5; 1; z^,^) + (v^tan((^^in/3) + l) F (1.5; 1.5; 2; 



vr^oCKn), (24) 



in which z^i^ = ^ — ^ cot and ({m,n) = m^, n^, i^^^y. for three cases in (0). Certainly, 
there is only half-width of stochastic layer, but the phase curves, laying out separatrix don't 
constitute an interest, being physically meaningless (they are unbounded on infinity). 



5 Summary and conclusions 

We have calculated the width of the stochastic layer around a separatrix, corresponding to a 
soliton solution of Korteweg-de Vries equation under hamiltonian perturbations. As far as the 
motion in this stochastic layer is chaotic ( in the meaning of mixing), so nonlinear wave solutions, 
which correspond to the phase curves within this stochastic layer will also posess certain chaotic 
properties. Let's now consider the question about the spatial region, where this chaotic behaviour 
can be registered. 

The time length, after which stochasticity in a nonlinear oscillator can be found is defined by 
Tc <C t, where Tc is a time of the decay of correlations. In the work |l^ the following estimation 
of Tc for a nonlinear oscillator is obtained: 

where i^' is a coefficient for the overlap of resonances, obtained in the previous section. As far 
as in our consideration time t for a nonlinear oscillator corresponds to the expression [x — ct) for 
nonlinear waves, so the region, where chaotic regimes for these waves can be detected, is: 

« (x - at) . (26) 
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Therefore, we can conclude that this region represents wave formation, which outstrips nonhnear 
wave or sohton, and propagates in the same direction with the same velocity. Under t = we 
obtain: 

(27) 

As far as with approach to a separatrix K —>■ oo, so for solitons the minimal distance on which 
chaotic behaviour should appear Xc —>■ 0, and this effect will be realized in a close vicinity of a 
soliton peak. The chaos, which we are searching for, will be manifested in the irregular small 
deflection from a smooth initial soliton profile. So, as a main result of this work, it can be inferred 
that long-period nonlinear waves and especially solitons of Korteweg-de Vries equation obtain 
chaotic properties in the presence of periodic hamiltonian external perturbations. 

Exactly this result was observed in an experiment with ion-acoustic and Langmuir waves 
in a nonmagnetized plasma There was registered a faint splash directly before the soliton. 
From the concepts developed in this paper this phenomenon can be explained in the following 
way. Together with soliton wave generator also produces a group of waves of small amplitude 
and almost zero frequencies. These waves are produced constantly with generator working, and 
therefore can be considered as a small deterministic periodic external perturbation. Stochastic 
destruction of soliton due to the influence of these waves, as it was described above, must realize 
itself just before the soliton peak, exactly what was registered in the experiment. 
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